Dynamical instabilities of dissipative solitons in nonlinear optical cavities with 

nonlocal materials 



OO 

o 
o 

(N 
> 

o 
O 
in 



> 

00 

o 

OO 

d 

OO 

o 



Lendert Gelens,^'El Damia Gomila,^'[l| Guy Van der Sande,^'^ Jan Danckaert,^ Pere Colet,^ and Manuel A. Matias^ 

^Department of Applied Physics and Photonics, Vrije Universiteit Brussel, Plemlaan 2, B-1050 Brussel, Belgium 
^IFISC, Instituto de Fisica Interdisciplinar y Sistemas Complejos (CSIC-UIB), 
Campus Universitat Illes Balears, E-07122 Palma de Mallorca, Spain 
(Dated: October 16, 2008) 

In this work we characterize the dynamical instabilities of localized structures exhibited by a 
recently introduced [Gelens et al, Phys. Rev. A 75, 063812 (2007)] generalization of the Lugiato- 
Lefever model to include a weakly nonlocal response of an intracavity metamaterial. A rich scenario, 
in which the localized structures exhibit different types of oscillatory instabilities, tristability, and 
excitability, including a regime of conditional excitability in which the system is bistable, is presented 
and discussed. Finally, it is shown that the scenario is organized by a pair of Takens-Bogdanov 
codimension-2 points. 
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I. INTRODUCTION 



Dissipative solitons (DS) are spatially localized struc- 
tures that appear in certain nonlinear dissipative media 
P, [3] . They have been found in systems such as chemical 
reactions 0, vegetation models 01, gas discharge sys- 
tems 0, fluids 0] and optics [1, i, [10, [H [H, [H, Q S ■ 
While in many instances these localized structures are 
stable, there are situations in which they develop differ- 
ent kinds of instabilities. Some instabilities lead to the 
formation of an extended pattern and therefore the local- 
ized character of the DS is destroyed. More interesting 
are the instabilities that, while preserving its localized 
character, induce the DS to start moving, breathing or 
oscillating [l^, [H, . Since DS can be considered 
as an entity on their own, these instabilities may lead 
to dynamical regimes that appear not to be present in 
the dynamical behavior of the extended system. In this 
context it has recently been reported that DS arising in 
a prototype model for optical cavities filled with a non- 
linear Kerr media may show excitable behavior, while 
locally the system is not excitable [13, [IJ. Thus, ex- 
citability can be an emergent property arising from the 
spatial dependence, which allows for the formation of lo- 
calized structures. In that situation excitability is medi- 
ated by a saddle-loop bifurcation and the whole scenario 
is organized by a Takens-Bogdanov (TB) codimension-2 
point. In parameter space the TB point is located in the 
asymptotic limit in which the model becomes equivalent 
to the Nonlinear Schrodinger Equation (NLSE). 

Since this excitability scenario is an emergent property 
of the spatial dependence of the system, it is particularly 
important to characterize how this scenario may change 
when the nature of spatial coupling is varied. In the 



Lugiato-Lefever model (22| considered in Refs. [20[, [21| 
the spatial coupling arises from optical diffraction in the 
paraxial approximation and is therefore accounted for by 
a Laplacian term. Here, we consider an extension of the 
model, including a mildly nonlocal term which extends 
the range of spatial interaction . This extension of the 
original Kerr model is suggested by the recent availabil- 
ity of metamaterials, allowing to design an optical Kerr 
cavity where layers of right- and left-materials are alter- 
nated. This provides the possibility to strongly decrease 
the diffraction strength in the resonator, such that higher 
order spatial effects (e.g. nonlocal effects) start to dom- 
inate the dynamical behavior of the DS. In this work, 
it will be shown that the additional spatial interaction 
term is able to shift the bifurcation lines such that now 
two Takens-Bogdanov points move from infinity to finite 
parameter values, acting as organizing centers of a richer 
dynamical behavior. 



II. MODEL 

We consider an optical cavity with a Kerr-type nonlin- 
earity, driven by a homogeneous, coherent optical light 
beam. This system was first introduced by Lugiato and 
Lefever to study pattern formation in a driven nonlin- 
ear, passive optical resonator [22]. In this work, we will 
study a more general equation, which includes a bilapla- 
cian term: 
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This model equation has been obtained in Ref. 
[2^ . [2^ . [25} to describe the temporal evolution of the 
slowly varying envelope of the electric field E{x,t) 
in a double-layered optical cavity. One layer of the 
cavity consists of a conventional right-handed mate- 
rial, while the other layer is an optical left-handed 
metamaterial. x — {x, y) represents the plane trans- 
verse to the propagation direction. Eq. ([T]) has been 
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obtained under the same approximations under which 
the Lugiato-Lefever equation is vahd, i.e., the slowly 
varying envelope approximation, weak nonlinearity, 
the paraxial limit and a nearly resonant cavity. The 
bilaplacian term appears when taking into account a 
linear and weakly nonlocal response of the left-handed 
metamaterial. Furthermore, it has been shown that in 
this double-layered cavity the ratio between a and f3 can 
be drastically altered by changing the relative lengths 
of both material layers. Recent advances in the field 
of metamaterials [la, [23|, potentially allow to find DS 
modeled by Eq. ([1]) in a wide area of the parameter space. 

The first term at the right-hand side of Eq. ([T|) models 
the cavity losses, 9 represents the cavity detuning with 
respect to the driving field Ei^, V\ = d'^/dx'^ + d'^/dy^ 
is the transverse Laplacian term due to diffraction, and 
the transverse bilaplacian term modeling the lin- 
ear, weakly nonlocal response. The cubic self-focusing 
nonlinearity is given by i\E\'^ E. Eq. ([T]) has a homo- 
geneous steady state solution Es = Eo/{l + i{9 — Is)), 
where Is = \Es\'^ [H, HI]- The homogeneous solution be- 
comes modulationally unstable at Ig — 1, which is also 
the case for the regular Lugiato-Lefever equation [2^ . 
but in Eq. (P) the instability arises with two characteris- 
tic wavenumbers 23] . From this modulational instability, 
a subcritical branch of DS appears. In the remainder of 
this manuscript, we use the background intensity Is, the 
detuning 9 and the coefficient of the bilaplacian term (3, 
that is a measure of the strength of the nonlocality, as 
our control parameters. Without loss of generality, we 
take a — 1. 
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FIG. 1: Phase diagram of DS in the nonlocal Kerr cavity for 
9 — 1.23. DS do not exist in Region I, below the saddle node 
line (SN). DS are stable between the saddle-node bifurcation 
(solid line - SN) and the Hopf bifurcation (dashed line - H), 
namely in Regions II, III, VI and VIII (see also Fig. [2} . In 
Regions IV and VII delimited by a Hopf and a saddle-loop 
bifurcation (dotted line - SL) there are stable oscillatory DS 
while the stationary DS are unstable. In Region V the static 
DS is unstable and the phase space generated after the de- 
struction of the limit cycle at the SL induces a regime of 
excitable DS. Region VI corresponds to a regime of condi- 
tional excitability with both stable and excitable DS. Finally 
in Region VIII one encounters tristability: a stationary and an 
oscillatory DS coexist with the homogeneous solution. Where 
the saddle-node bifurcation line and the Hopf bifurcation line 
meet is a Takens-Bogdanov (TB) bifurcation. 



III. DYNAMICAL BEHAVIOR 
A. Preliminary remarks 

In this section, we show the different possibilities of 
dynamical behavior of DS in parameter space (see Fig. 
[T]). Since DS are radially symmetric they correspond to 
stationary solutions of the radial form of Eq. ((T]) with 
boundary conditions drE{r = 0) = and drE{r — > 
oo) = 0. This equation can be solved numerically us- 
ing a Newton method [2l|, IH, |2^ . This approach is very 
accurate and automatically generates the Jacobian op- 
erator, whose eigenvalues determine the stability of the 
solutions. Note that this method finds both stable and 
unstable stationary solutions. DS can undergo two kinds 
of instabilities, radial instabilities which preserve the lo- 
calized character of the structure and azimuthal instabil- 
ities which lead to the formation of extended patterns. 
The last ones appear only for large values of the back- 
ground intensity (/« close to 1). We focus here on the 
radial instabilities, so phase diagrams (Figs. [T] and [2) are 
plotted only up to Is — 0.93, before the azimuthal insta- 
bilities take place. 

Since the system has three parameters (Is, 9, (3), for 



the sake of clarity we fix the detuning at 6* = 1.23 in 
this section and analyze a slice of the whole parameter 
space. Fig. [1] shows the region of the parameter plane 
{P, Is) that contains the most relevant regimes of dynam- 
ical behavior of the system. The line that dominates this 
parameter plane has the shape of a deformed parabola, 
and is a line of saddle-node (SN) bifurcations in which 
two DS are created. Below this line one has Region I 
where no DS exist. We recall that in all the parame- 
ter range covered by Fig. [T] the spatially homogeneous 
solution is always stable. The different regimes above 
this line are organized by two codimension-2 Takens- 
Bogdanov points, TBi (/3 = 0.00987, Is = 0.7741) and 
TB2 (/3 = 0.02944,/,. = 0.6707), discussed in the next 
two Sections. 

Without the bilaplacian term there is only one TB 
point [l^jUlj. In that case, the TB point is found only 
asymptotically for the limit ^ ^ 00. As wc will discuss 
later, that TB point corresponds in fact to the TBi point 
found here. Therefore, the nonlocality has brought this 
bifurcation to finite parameter values allowing us to fully 
study the different dynamical regimes around the TBi 
point. 
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FIG. 2: Zoom of Fig.[T]near TBi (a) and TB2 (b). Lines have 
been determined using the method explained in Section Till Al 
with the exception of the SL hne in panel (b). In that case 
the filled squares display the location of the SL obtained from 
numerical integration of Eq. ([T}, while the grey line through 
these points is only to guide the eye. 9 = 1.23. 



B. Dynamical scenario around the TBi point 

A Takens-Bogdanov (or double-zero) bifurcation is as- 
sociated to the presence of two (non-diagonalizable) de- 
generate null eigenvalues [H, [13] ■ Such a bifurcation oc- 
curs when, in a line of SN bifurcations, one of the modes 
transverse to the center manifold (of the SN bifurcation) 
passes through zero, implying that this transverse mode 
switches from stable to unstable or vice versa. If this 
transverse mode is stable we will denote the SN bifurca- 
tion line as SN~ , while we use SN+ if this mode is unsta- 
ble. Throughout the remainder of this paper, we will use 
~ for bifurcations for which there is a stable emerging so- 
lution and + if the emerging solutions are unstable. H~ 
will describe a supercritical Hopf, and a subcritical 
one. 

Another feature of a TB point is that two new bifur- 
cation lines emerge from it [31]: a Hopf bifurcation line 
[36(] and a saddle- loop (homoclinic) bifurcation line [37[. 
In order to specify whether the cycle that emerges from 
the saddle-loop bifurcation is stable or unstable, it is use- 
ful to define the saddle-quantity v. For low-dimensional 
dynamical systems this quantity is given by = As + A„, 
with Au > and A^ < the unstable and stable eigenval- 
ues of the saddle, respectively. The emerging cycle will 
be stable if z^ < 0, and unstable in the opposite case: 

> [H. 




FIG. 3: Qualitative evolution of the phase space at la ~ 0.8 
and 9 = 1.23 for increasing values of P, corresponding to the 
different dynamical behavior in a horizontal cut of Fig. [T] 
From (a) to (j), one goes from stable DS (Region II) to an 
oscillating DS (Region IV), followed by an excitable DS (Re- 
gion V), a conditional excitable DS (Region VI), a coexistence 
of stable DS and oscillatory DS (Region VIII) and finally a 
stable DS (Region III). 



At the TBi point, the saddle-node bifurcation is stable 
(SN~), the Hopf is supercritical (H^), and the saddle- 
loop creates a stable cycle (SL^). Along the SN^ line 
a pair of stationary DS are created, one stable (upper 
branch) and the other (middle branch) unstable along a 
single direction in phase space (thus, a saddle point in dy- 
namical systems parlance). So, Region II is characterized 
by stable DS coexisting with the spatially homogeneous 
solution. A qualitative sketch of the most relevant differ- 
ent kinds of behavior in the system as /3 is increased can 
be found in Fig. [3] (this corresponds to a horizontal line 
in Fig. [2] for /g « 0.8). In this figure, panel (a) reflects 
the behavior inside Region II, where the DS is the stable 
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focus, the homogeneous solution is the stable node and 
the middle branch DS is the saddle. 

The upper branch DS solution becomes unstable at the 
supercritical Hopf bifurcation, H~, and leads to Region 
IV, characterized by oscillatory DS, i.e., autonomous os- 
cillons. Fig. [3] (b) illustrates the behavior past the Hopf 
bifurcation: the stable oscillatory DS and the unstable 
focus in the center can be seen. Approaching the SL~ 
line a saddle-loop (homoclinic) bifurcation takes place, 
in which the limit cycle (oscillatory DS) becomes a ho- 
moclinic orbit of the saddle (middle branch DS). The SL^ 
is a global bifurcation, and cannot be detected through 
a local analysis. Thus, in this study it has been deter- 
mined through numerical simulations of Eq. ([T]). Panel 
(c) in Fig. [3] illustrates the cycle growing in amplitude 
and approaching the saddle, while in panel (d) the cycle 
has become the homoclinic orbit of the saddle. The ap- 
proach of the stable cycle to the saddle can also be seen 
quantitatively in panel (a) of Fig. 3] in this figure bifur- 
cation diagrams corresponding to vertical cuts in param- 
eter space, i.e. with /3 fixed (cf. Fig. [J), are presented. 
Beyond the SL~, the behavior of the system is excitable 
(Region V) [13, Hlj , in particular of type (or class) I [s^l , 
as the excitability threshold is the stable manifold of the 
saddle. An excitable excursion is achieved when local- 
ized perturbations beyond this threshold are applied to 
the spatially homogeneous solution. Fig [3] (e) sketches 
the phase space in the excitable regime past the saddle- 
loop bifurcation. 

Inside Region V, no stable DS are found, only unsta- 
ble DS solutions exist, and the system is excitable for all 
values of Is above the SN"*" line. This line which is just 
the continuation of SN^ past the TBi point was is not 
observed in [l^, HH since the TB point without bilapla- 
cian term is located at — > oo. The SN+ line creates a 
saddle-unstable node pair (middle and upper branch, re- 
spectively). The pair of unstable solutions can be seen in 
the bifurcation diagram shown in Fig. correspond- 
ing to a vertical cut of the parameter plane just to the 
right of TBi point. 



C. Dynamical scenario around the TB2 point 
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An even richer scenario is found around the TB2 point. 
Comparing Fig. [2Ja) and (b) , one can see that the TB2 
does not yield the same scenario as around TBi. In TBi 
the two lines that emerge involve stable objects (II~ and 
SL~ create and destroy, respectively, a stable limit cycle), 
while in TB2 the Hopf line is subcritical (involving an 
unstable cycle). Furthermore, in the TBi point the two 
lines, H~ and SL~ are tangent to the line SN~ of stable 
saddle-node bifurcations while the opposite happens for 
TB2 (the lines unfolding are tangent to SN"*"). 

These differences correspond to a change of sign in a 
term of the normal form of the TB bifurcation 13 ill : 
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FIG. 4: Bifurcation diagrams showing the maximum intensity 
of the DS as function of the background intensity Is , for: (a) 
13 = 0.005; (b) /3 = 0.015; (c) (3 = 0.0268; (d) /? = 0.0272; (e) 
P = 0.05. The lowest solid line represents the stable homoge- 
neous solution. The lowest dashed line shows the maximum 
intensity of the unstable middle branch DS. Above that line, 
the upper branch DS is shown, where the solid line stands 
for the stable DS and the dashed line for the unstable DS. 
In (a), (c) and (d) the system exhibits stable oscillatory be- 
havior (see text). In these cases, the maximum and minimum 
intensity of the oscillating DS are depicted as filled squares. 
The grey dashed line, representing the unstable limit cycle, is 
only a guiding line for the eye. When the limit cycle touches 
the middle branch DS, a saddle-loop bifurcation occurs. 9 = 
1.23. 
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The case that is most often discussed in the hterature 
leads to the scenario discussed in Section fill B[ and cor- 
responds to s = — 1. A supercritical Hopf line unfolds 
from the TB point tangent to the SN~ segment, and a 
saddle-loop, in which a stable cycle is destroyed unfolds 
tangent to the Hopf line. In the s = -1-1 case, however, 
a subcritical Hopf bifurcation line that creates an unsta- 
ble cycle unfolds tangent to the SN+ segment, and the 
unstable cycle is destroyed in a saddle-loop also tangent 
to the Hopf bifurcation. This can be seen noticing that 
changing the sign of s is equivalent to perform the sub- 
stitutions t i-^ ~t and y i— > — This scenario fits nicely 
with what is observed around the TB2 point, and, thus, 
TBi and TB2 correspond to the two possible cases of Eq. 
(I2|), with s = — 1 and s = -fl, respectively. 

Crossing the H+ line coming from Region V, the unsta- 
ble (upper branch) DS exhibits a subcritical Hopf bifur- 
cation and becomes a stable focus. The cycle that is cre- 
ated is unstable and has only a single unstable direction. 
In Region VI the system is then bistable (upper branch 
DS and homogeneous solutions coexist) and the unstable 
cycle is the basin boundary of the upper branch stable 
DS as qualitatively illustrated in Fig. [3] (f). Just after 
the bifurcation this basin of attraction is small (since the 
initial cycle amplitude is zero) and grows as one moves 
away from the bifurcation line. This can be seen in the 
bifurcation diagram displayed in Fig. [DJc) corresponding 
to a vertical cut of Figs. [T]and[2]at /3 — 0.0268. The upper 
branch DS becomes stable around Is = 0.8. The dashed 
grey line has been drawn with the purpose of guiding 
the eye only. It represents the unstable cycle that we do 
not compute. This bifurcation diagram also shows the 
existence of a stable limit cycle (plotted as squares) cor- 
responding to an oscillatory DS. The stable limit cycle 
comes from a fold of cycles (FC) bifurcation, discussed 
in more detail later, which takes place at a larger value 
of Is- Decreasing Is the stable limit cycle disappears 
at a saddle-loop bifurcation which takes place when the 
stable limit cycle becomes the homoclinic orbit of the 
saddle (middle branch soliton). In Fig. \Ms) Region VI 
corresponds to the values of Is limited on the left by the 
subcritical Hopf (H+) and the right by the saddle loop 
(SL+ and SL~). Precise values for the SL~ obtained 
from numerical integration of Eq. ([T]) are plotted as filled 
squares in Fig. [2tb) while the grey SL line joining the 
points have been drawn to guide the eye. As we will dis- 
cuss in the next subsection, Region VI corresponds to a 
regime of conditional excitability. 

In Fig. HKc), the FC on the one hand and the SL^ 
on the other limit a new region of tristability where a 
stationary DS, a oscillatory DS and the homogeneous 
solution coexist. In Fig. [2Kb) the tristable region is la- 
beled as VIII. In Fig. [3] the panel (g) sketches the phase 
space at the saddle-loop while the panel (h) illustrates the 
tristable regime. Increasing /3 or Is in parameter space 
the stable cycle decreases in amplitude while the unstable 
limit cycle increases until both the stable and unstable 
cycles are destroyed in the fold bifurcation (Fig. EJi)). 
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FIG. 5: Dynamical evolution of the system in the tristable 
regime, where the spatially homogeneous fundamental branch 
(not shown), stable oscillatory DS regime (solid line), and 
stable focus (gray line) coexists. 




FIG. 6: A qualitative sketch of the phase space at Is ~ 0.9 
and 6 — 1.23 for increasing values of /3, corresponding to a 
horizontal cut of Fig.[2Ib), is depicted in Fig. |S] Increasing P 
from (a) to (d), one starts with excitable DS (Region V), going 
through a region of bistability between a stable limit cycle and 
the stable homogeneous solution (Region VII), followed by a 
coexistence of stable DS and oscillatory DS (Region VIII) and 
finally a stable DS (Region III). 

Figure [5] shows the time evolution in the tristable regime 
obtained starting from an initial condition belonging to 
the basin of attraction of the limit cycle and from an ini- 
tial condition within the basin of attraction of the stable 
DS. 

In the phase diagram (Fig. [^Jb)) the SL^ line can also 
be located to the left of H+ (Region VII) . The bifurcation 
diagram in this case would be similar to the one shown 
in Fig. UKc), except for the fact that the SL~ becomes a 
SL+ and the H"*" takes place at a larger value of Is than 
the SL~. Therefore VII is a region of bistability where a 
stable limit cycle corresponding to the oscillatory upper 
branch DS coexist with the homogeneous solution while 
the steady state upper branch DS is an unstable focus. 
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A qualitative sketch of the phase space when going from 
region V through Vll and Vlll is depicted in Fig. [HI 



The last region in the phase diagram to be described 
is Region III, located to the right of the SL+ and FC 
lines. The phase space corresponding to this broad region 
is sketched in Figs, and [DJd). The upper branch 
stationary DS is a stable point and coexists with the 
homogeneous solution. Fig. IDJe) shows a quantitative 
bifurcation diagram corresponding to a vertical cut of 
Fig.[I]for P = 0.05, to the right of the TB2. For < 0.7 
only the homogeneous solution exists (Region 1) while 
above Is —0.7 one enters in Region 111 of coexistence of 
the stable DS with the homogeneous solution. 



We now analyze in detail the fold of limit cycles bifur- 
cation. It appears at a secondary codimension-2 point 
named resonant side-switching, where the fold and the 
saddle-loop bifurcation, that delimit the region of trista- 
bility, coalesce [H, [s^l- This occurs, roughly, around 
P = 0.0272 and = 0.82. Fig. H^d) depicts a bifurca- 
tion diagram close to this codimension-2 point, for which 
the region of existence of stable limit cycles is quite nar- 
row: the represented square is very close both to the 
saddle-loop bifurcation and to the fold of cycles. In the 
phase diagram shown in Fig. this codimension-2 cor- 
responds to the point where the FC and the SL (SL~ and 
SL+) lines meet. 
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FIG. 7: The spectrum of the middle-branch DS at the saddle- 
loop bifurcation for (a) /? = 0.0268 and Is = 0.84, and (b) 
f] = 0.0294 and 1^ = 0.671. The filled black circles are the 
eigenvalues corresponding to the modes with a localized spa- 
tial profile, e = 1.23. 



The occurrence of the resonant side-switching bifurca- 
tion is related to the eigenvalue spectrum of the saddle, 
namely, the saddle quantity of the middle branch DS ap- 
proaches zero, and the cycle emerging from the saddle- 
loop changes from " stable" to " unstable" . Close to TB2 
point the saddle-loop bifurcation must destroy an unsta- 
ble cycle (transition from Region VI to 111), hence > 0, 
while after the fold has taken place it destroys a sta- 
ble cycle (transition from Region Vlll to VI), so < 0. 
Fig. [7] shows the spectrum of the middle branch DS for 
two different sets of parameters, one after the formation 
of the fold (saddle-loop bifurcation with a stable cycle, 
top panel) and one very close to the TB2 point (saddle- 
loop bifurcation with an unstable cycle, bottom panel). 
The eigenvalues of the modes relevant for the dynamics 
are highlighted with filled circles, and indeed, the saddle 
quantity is positive > 0, albeit quite small {i^ ~ 0.003) 
in the bottom panel, compatible with the fact that it in- 
volves an unstable cycle. One can observe that close to 
the SL+ that emerges from the TB2 point, where the sad- 
dle quantity is positive, four localized modes play a role 
in the dynamical behavior of the system. However, when 
moving away from the SL"*" line, only two localized modes 
determining the dynamics remain. So, in conclusion, we 
can say that the system is essentially two-dimensional for 
the same reasons as discussed in Rcf. 2l], except close 
to the SL"^ line, originating from the TB2 point. 



D. Excitability and conditional excitability 

When crossing a saddle-loop bifurcation from a region 
where a stable oscillatory DS exist one enters in a ex- 
citable regime. As discussed in Subsection IIII Bl this 
scenario takes place when going from Region IV (where 
the oscillatory DS is stable) to Region V where DS ex- 
hibit excitability. Of course it is also possible to enter in 
Region V from the other side, namely from the side of 
the TB2 point. In fact a similar scenario is found when 
going from Region VII where the oscillatory DS is sta- 
ble to Region V. In any case the stable manifold of the 
saddle (middle branch DS) plays the role of excitability 
threshold, so the excitable response is triggered only by 
localized perturbations of the homogeneous solution that 
bring thesystem beyond this threshold. Excitability is of 
Class I [32] , characterized by long response times for per- 
turbations that leave the trajectory close to the saddle 
in phase space. 

Close to the TB2 point there is another region, VI, 
where one can enter crossing a saddle-loop line. How- 
ever, the dynamical behavior in Region VI is qualitatively 
different from Region V since the upper branch LS is an 
unstable focus in the last one, while it is stable in the for- 
mer. As discussed before, the upper branch LS has been 
made stable by the subcritical Hopf bifurcation H+ that 
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separates Region V from VI. The H"*" bifurcation gen- 
erates also an unstable limit cycle, which is not present 
in Region V. Therefore, although the transition from Re- 
gions III or VIII to Region VI goes through a saddle-loop 
bifurcation, the scenario must be qualitatively different 
from the one discussed above which leads to the usual 
excitability found in Region V. In Region VI, one finds 
a regime of conditional excitability, in which the DS is 
simultaneously excitable and bistable. The excitable be- 
havior is also Class I. 

To clarify what conditional excitability means, we refer 
to the phase space sketched in Fig.[3]Jf). In this situation, 
while as usual, perturbations of the homogeneous solu- 
tion that are not able to cross the excitability threshold 
(stable manifold of the saddle) lead to normal relaxation, 
there are two possible different dynamical responses for 
supra-threshold perturbations. If a localized perturba- 
tion of the homogeneous state brings the system inside 
the basin of attraction of the stable focus, namely inside 
the unstable cycle, the system jumps from the fundamen- 
tal solution to this attractor. Therefore after this pertur- 
bation the system relaxes to the stable DS in a oscillatory 
way. The grey line in Fig. [5] shows the dynamical evolu- 
tion of the maxima of the peak in this situation. Instead, 
for localized perturbations of the homogeneous solution 
which bring the system beyond the stable manifold of 
the saddle but outside the unstable cycle, the response 
is excitable. The system exhibits a large response corre- 
sponding to a circulation around the unstable limit cycle 
before returning to the stable homogeneous solution. The 
black line in Fig. [S] shows the time evolution of maxima 
of the peak for a excitable trajectory. 

So, in summary, the dynamical response of pertur- 
bations is more complex than simply sub- and supra- 
threshold, and for the latter type of perturbations two 
possible regimes are possible. When going from Region 
VI to Region V at the 11+ line, the upper branch DS be- 
comes unstable and the unstable limit cycle responsible 
for the conditional excitable response to supra-threshold 
perturbations disappears, so the conditional excitability 
becomes a usual one. 
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FIG. 8: The conditional excitability regime is illustrated for 
P = 0.0272 and /« = 0.8. The solid line (excitable trajectory) 
corresponds to time evolution of the maxima of the peak after 
a localized perturbation of the homogeneous state that brings 
the system beyond the excitability threshold but outside the 
basin of attraction of the upper branch DS. The grey line 
corresponds to a stronger localized perturbation that brings 
the system inside the boundary of attraction of the upper 
branch DS. 
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FIG. 9: Location {9tb ,I3tb) of the two TB points as the 
detuning is increased. The solid line and the dashed line cor- 
respond to the TBi and TB2 point, respectively. 



IV. CONCLUDING REMARKS 

We have studied the nonlinear dynamical behavior of 
2D localized structures in a model for an optical cav- 
ity filled by a Kerr nonlinear medium and a left-handed 
metamaterial [23|. The model is a generalization of the 
Lugiato-Lefever equation 22], and includes higher order 
spatial effects arising from the weakly nonlocal response 
of the metamaterial. In this system, we have shown the 
existence of regions with stationary, oscillating and ex- 
citable localized structures. Furthermore, we have shown 
that the different bifurcation lines originate from two 
Takens-Bogdanov (TB) codimension-2 points, which is 
a strong signature for the presence of a homoclinic bi- 
furcation. This homoclinic bifurcation offers a route to 



excitable behavior of the 2D localized structures. Finally 
an extra secondary codimension-2 point (resonant side- 
switching bifurcation) creates a fold of cycles that leads 
to two new regimes, one of tristability and one of condi- 
tional excitability. 

Without the nonlocal terms the different dynamical 
regimes of the DS were organized by a TB point lo- 
cated at the limit of infinite detuning, where the Lugiato- 
Lefever equation reduces to the (conservative) Nonlinear 
Schrodinger Equation (NLSE) 20, 21j. Here, we demon- 
strate that the presence of higher order spatial effects 
brings two TB points to finite parameter values. In Fig. 
[51 we provide further evidence that TBi is unfolding from 
the NLSE at the hmit (3^0 and then 6* 00 of Eq. 
([I]). This hypothesis is supported by the fact that 2D 
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solitons in the NLSE have at least a twofold degeneracy 
[sH - Furthermore, Fig. [5] also gives some evidence that 
TB2 comes from a certain conservative limit at 6* ^ cx), 
/3 ^ of Eq. However, if this is a singular limit, this 
does not necessarily imply that TBi and TB2 have to 
evolve towards the same point in the conservative limit. 
The numerical evidence given in Fig. [9] is of course not 
fully conclusive, and this unfolding of both of the TB 
points is presently under investigation. 
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